Cu Nuclear Quadrupole Resonance Study of the Spin-Peierls Compound 
Cui-^Mg^GeOa: A Possibility of Precursory Dimerization 
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We report on a zero- field ^''Cu nuclear quadrupole resonance (NQR) study of nonmagnetic Mg impurity 
substituted Cui-^MgajGeOs (single crystals; the spin-Peierls transition temperature rsp~14, 13.5, and 11 
K for a:=0, 0.0043, and 0.020) in a temperature range from 4.2 K to 250 K. We found that below r*-77 
K, Cu NQR spectra are broadened and nonexponential Cu nuclear spin-lattice relaxation increases for 
undoped and more remarkably for Mg-doped samples. The results indicate that random lattice distortion and 
impurity-induced spins appear below T*, which we associate with a precursor of the spin-Peierls transition. 
Conventional magnetic critical slowing down does not appear down to 4.2 K below Tsp. 
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The discovery of the first inorganic spin-Peierls com- 
pound CuGeOs (the transition temperature Tsp^ld 
K) [y and subsequent reports on unprecedented impu- 
rity effects [^-0 have renewed the interests of a quasi- 
one-dimensional spin 5=1/2 Heisenberg antiferromagnet 
coupled to phonons. No soft mode of phonon at T^p is 
one of the characteristics of CuGeOa. An appreciable 
inter-chain exchange interaction g], lattice and phonon 
anomaly perpendicular to the chain [B 10| are different 
from an ordinary spin-Peierls system or conventional the- 
oretical result |]11|. A pseudogap in the magnetic excita- 
tion spectrum below 20 K observed by inelastic neutron 
scattering |1J] and a local dimerization until at least 40 
K observed by diffusive X-ray scattering ^M resemble a 
pseudogap of the electronic Peierls materials |l4| . In con- 
trast to conventional competition between Neel ordering 
and lattice dimerization, the impurity substitution (Zn, 
Mg, Ni, or Si) for Cu or Ge induces a dimerized antifer- 
romagntic ordering state P-H|, where a spin- wave mode 
coexists in the spin-Peierls gap [Q. The coexistence at 
r=0 K is understood within the framework of the phase 
Hamiltoninan |]l5| . 

Nuclear quadrupole resonance (NQR) and nuclear 
magnetic resonance (NMR) are unique and powerful 
techniques to study low frequency dynamics and local 
spin fluctuations in space. The intensive studies using 
Cu NQR and NMR techniques have revealed many as- 
pects of CuGeOa |l6|-[l9||. The spin gap opening at Tgp 
was evidenced by an abrupt decrease of the Cu nuclear 
spin-lattice relaxation rate l/Ti without any appreciable 
change of Cu NQR spectrum ||l^, although a singlet- 
triplet excitation, a spin gap and ion displacement were 
directly confirmed by neutron scattering ||2^,|2^. The 
above Tsp spin dynamics, the Cu l/Ti, is understood by 
a quasi-one-dimensional 5*= 1/2 antiferromagnctic corre- 
lation without spin-phonon coupling [n6U23j23]. To our 



knowledge, however, there are no reports of Cu NQR 
spectrum far above Tsp for CuGeOa nor of impurity ef- 
fects on the low frequency spin dynamics, after the work 
on Zn doping |2|] . 

In this Letter, we report the high temperature mea- 
surement of Cu NQR spectrum for undoped CuGeOa 
and the Cu NQR study of nonmagnetic Mg impurity 
substitution effect on Cui_,2,Mgj,Ge03 (single crystals; 
Tsp ~14, 13.5, and 11 K for a;=0, 0.0043, and 0.020) in a 
wide temperature range of r=4. 2-250 K. We found that 
the broadening of Cu NQR spectra and the increase of 
nonexponential Cu nuclear spin-lattice relaxation occur 
below about 77 K much higher than Tsp, which suggest 
a precursor of the spin-Peierls transition. We did not ob- 
serve critical divergence of l/Ti for x=0.020 down to 4.2 
K, although the magnetic ordering occurs at about 2.5 
K|. ^ 

The single crystals grown by a floating-zone method 
are well characterized in Ref. [0. Zero-field ^^Cu 
NQR spin-echo measurements were carried out with 
a coherent-type pulsed spectrometer. NQR frequency 
spectra with quadrature detection were measured by in- 
tegration of the ^^Cu nuclear spin-echoes as the frequency 
was changed point by point. Nuclear spin-lattice re- 
laxation was measured by an inversion recovery spin- 
echo technique, where the ^'^Cu nuclear spin-echo am- 
plitude M{t) was recorded as a function of time in- 
terval t, between an inversion vr-pulse and a 7r/2-pulse 
(tt — t — 7r/2 — TT-echo). 

Figure 1 shows ^'^Cu NQR spectra for undoped x—Q 
(a) and for Mg doping of a;=0.020 (b) in the temper- 
ature range of r=4.2-250 K. The observed Cu NQR 
spectra for Mg doping are nearly symmetrically broad- 
ened, not of Gaussian nor of Lorentzian type but rather 
have a triangle-shaped line profile for a;=0.020 at 4.2 
K. Implication of the characteristic line shape is not 
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FIG. 1. The temperature dependence of ^^Cu NQR fre- 
quency spectrum for x—0 (a) and for a;=0.020 (b). The tem- 
perature dependence of the peak frequency ^^v (c) and of the 
hnewidth (fuU-width at half maximum) ^^Ai/1/2 (d) for x=0, 
0.0043. and 0.020. The dashed hues in (d) are guides to the 
eye. 



clear. In general, the Cu NQR frequency v is given 
by J^ = (e^g(5/2/i)-\/l + ri'^/3, where eq is the maximum 
component of the electric field gradient tensor at the nu- 
clear site, Q is the nuclear quadrupole moment, and 77 
(0 < 7? < 1) is an asymmetry factor |25[ . The muon spin 
relaxation measurements have not detected any static in- 
ternal magnetic field for Zn-doped samples above about 
4 K 0. Thus, it is likely that the random distribution 
of the electric field gradient (eg and rj) is the origin of 
broadening. 

Figure 1(c) shows the temperature dependence of the 
peak frequency ^^v for x^O, 0.0043, and 0.020. The 



observed linear temperature dependence of 



63, 



nearly 



independent of Mg content, is similar to that of the 
cuprate mono-oxide CuO [2^. Figure 1(d) shows the 
temperature dependence of the linewidth defined as full- 
width at half maximum ^'^Ai^i/2. The linewidth as well 
as the peak frequency does not show any appreciable 
change at T^p, in agreement with the previous reports on 
CuGeOa |l|@|2|l. The linewidth of the x=0.020 sample 
is about 3 times larger than that of x—0 at T > T* . In 
Fig. 1(d), we obtain a significant result that ^'^Azyi/2 in- 
creases rapidly below about 77 K >> Tgp (denoted as T*) 
for all samples including CuGeOa. For each x, ^•^Ai/i/2 
at 4.2 K is about 2 times larger than that above T*. 

NQR is a measure of deviation of charge distribution 
from cubic symmetry around the nuclear site, being quite 
sensitive to crystal imperfections. The observed nearly 
symmetric line shape implies a random distribution of 
local charge pq]. If the origin of ^^Ai/i/2 is a static distri- 
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FIG. 2. Mg-doping dependence of ^''Cu recovery curves 
1-M(i)/M(oo) at 10 K (a) and at 16 K (b) for x=0, 0.0043, 
and 0.020. The solid curves are the least-squares fitting re- 
sults using eq. (1). 



bution of lattice distortion around crystal imperfection, 
^^Ai/i/2 would decrease as the temperature is decreased 
so as to scale with the temperature dependence of ^^v. 
However, the actual ^■^Ai/i/2 increases below T*. Thus, 
the inhomogeneity of lattice distortion must depend on 
temperature and must increase rapidly below T*. The 
pretransitional lattice fiuctuations above Tsp observed by 
the diffraction experiment P4 I, which are explained by 
the random-phase approximation calculation pq |, may 
be closely related with the increase of ^■^Aj/i/2. Accord- 
ing to the recent quantum Monte Carlo simulation p9[ , 
a precursory dimerization takes place near the edges far 
above T^p. T* corresponds to the onset of preformed 



L Sp 

dimer bonds 



In the soliton picture ]30||, Tsp is an order- 



to-disorder transition temperature of locally dimerized 
segments. T* may correspond to the onset of develop- 
ment of the interchain correlation between the soliton 
and antisoliton. 

Figure 2 shows Mg-doping effect on ^'^Cu recovery 
curve p{t) = 1 — M{t)/M{oo) of the ^'^Cu nuclear mag- 
netization M{t) at 10 K (< Tsp) (a) and at 16 K (> Tsp) 
(b). The recovery curve changes from a single exponen- 
tial function to nonexponential one as Mg is substituted 
and as the temperature is decreased. To account for the 
nonexponential function, we assume a minimal model, 
which consists of a host homogeneous relaxation process 
and of a single inhomogeneous relaxation one. The solid 
curves are the least-squares fitting results using the fol- 
lowing equation [ pi[ . 

Pit) = p(0)exp[-(i/Ti)ArQfl - ^tJT,]. (1) 
The fit parameters arep(O), (Ti)nqr @ and ti. p{0) 
is a fraction of an initially inverted magnetization, and 
{Ti)nqr is the nuclear spin-lattice relaxation time due to 
the host Cu spin fluctuations, ri is an impurity-induced 
nuclear spin-lattice relaxation time, which is originally 
termed a longitudinal direct dipole relaxation time, be- 
cause the second term of eq. (1) is derived from a ran- 
dom Ti process of 1/Ti(r) — C/r^ (C is a constant, and 
r is a distance between an impurity-induced spin S and 
a Cu nuclear spin /) through a direct dipole coupling 
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FIG. 3. Log-log plots of ^^(l/ri)jvQ_R (a) and ''''(l/n) 
(b) as functions of temperature for undoped x=0 and for 
Mg-doped samples of x=0.0043 and 0.020. In (a), the solid 
line is a T^ function, and the dashed curve is the least-squares 
fitting result using a function of l/Ti = i?iexp(— 2A/r). 



(xI±Sz/r^ {Sz is the z-component of S, and I± is a rais- 
ing or lowering operator of /) . The randomly distributed 
impurity-induced spins yield the stretched exponential 
function of eq. (1). The original Mg ion does not carry 
spin 1/2. Atomic defects or Mg ions cut chains into seg- 
ments. The existence of spatially extended staggered mo- 
ment induced by an edge or an impurity has been pointed 
out for a finite or a semi- infinite chain P9| , p0| , p3[ . The 
assumption of impurity-induced spins could be only a 
working hypothesis to introduce ti. Since the essence of 
the stretched exponential function is randomness in the 
Ti process, one may speculate that Ti{r) with a local 
spin density induced by Mg is approximated by a power 
law, leading to the stretched exponential function. In the 
soliton picture l^^l ? the soliton which stays in the middle 
of a segment or near the edges due to an interchain cou- 
pling, carries spin 1/2, so that it can act as an impurity 
spin. 

Figure 3 shows log-log plots of ^^{l/Ti)]yQji (a) and 
^^(1/ti) (b) as functions of temperature for x—0, 0.0043 
and 0.020. Far above T*, '^^{1/Ti)nqr for Mg-doped 
samples is nearly the same as that for an undoped one. 
For undoped and x=0.0043, the activation-type tem- 
perature dependence of ^'^(l/ri)jvQ_R is observed below 
Tsp. The spin gap A is estimated to be ~26 K by fits 
of l/ri=i?iexp(-2A/r) {Ri and A are fitting parame- 
ters) |3J], which agrees with the value estimated from 
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FIG. 4. The temperature and Mg-doping dependence of 
^''{1/TiT)nqr- The solid curve is the squared static uniform 
susceptibility xt, of undoped CuGeOa reproduced from Ref. 
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the static susceptibility ^. For a;=0.020, however, the 
temperature dependence of ^^{1/Ti)nqr is changed into 
a power-law type ('^T^), probably because of an inhomo- 
geneous distribution of A(r). Conventional critical diver- 
gence toward the magnetic ordering does not appear. 

Below around T*, ^'^(l/ri) immediately increases as 
the temperature is decreased down to Tgp even for a::=0, 
which indicates the increase of the impurity-induced spin 
correlation. Far below Tsp, ^^(1/ri) is systematically en- 
hanced by Mg doping, which is due to an increase of the 
number of impurity relaxation centers. It is likely that 
the origin of ^^(1/ri) for a;=0 is due to non-intentionally 
introduced imperfections (defects, dislocations, . . . ). The 
actual sample is not a perfect crystal, because the ob- 
served linewidths of ^'^Cu NQR spectra of our x—0 (^180 
kHz at 4.2 K, ~100 kHz above 100 K) are broader than 
those expected from Ti or T2 broadening (a few kHz), i.e. 
inhomogeneous broadening. The estimated ^^At^i/2 for 
x—0 is nearly the same as or somewhat sharper than the 
reported values below 40 K [16, 17, 22, 27]. In Fig. 3(b), 
^^Ti for each x makes a kink at around Tgp. In terms of 
an impurity spin picture, ^'^(1/ti) is nearly proportional 
to the life time of the impurity spin scattered by the host 
magnetic excitations. Then, the kink of ^^(1/ri) reflects 
a change in the host magnetic excitation spectrum at 
the true transition temperature Tsp, which is evident in 

''HI/Ti)nqr. 

Figure 4 shows the Mg-doping effect on ''^{l/TiT)NQn. 
For comparison, the squared static susceptibility Xu of 
undoped CuGeOa (solid curve) is also reproduced from 
Ref. [1]. In general, 1/TiT is the low frequency dy- 
namical spin susceptibility at an NQR frequency summed 
over a momentum space via a nuclear-electron coupling 
[35]. For a;=0 above Tsp, '^^{1/TiT)nqr is understood 
by the sum of the staggered spin susceptibility x{q = 



7r)~l/T and the Bonner-Fisher- type uniform spin sus- 
ceptibiUty Xu (to be exact, x«) lUiIll- The upturn of 
^^{1/TiT)nqr just above Tgp, which is ascribed to the 
staggered xil = ''r)~l/T, is suppressed by Mg doping of 
x=0.020. Then, the contribution from the uniform mode 
q^O is uncovered, being similar to the temperature de- 
pendence of Xu- The q = TT mode is easily affected by 
imperfection, comparatively more than the q—0 mode, 
as can be seen for La2Cui_2:Zna;04 [|6|. However, one 
should note that ^^il/TiT)NQR for x=0.020 and xl for 
x=0 above T^p are similar but do not completely agree 
with each other. The suppression of ^^(l/Tir)jvQ_R for 
a;=0.020 begins from below 60-120 K more steeply than 
that of Xu below about 50 K. Thus, the further mecha- 
nism of the suppression is needed. The precursory dimer- 
ization enhanced by Mg below around T* is a possible 
candidate. Our observations of the suppressed xi<l — ^) 
and of the deviation between ^^{l/TiT)]yQji and Xu foi' 
Mg doped samples will be constraints on dynamical the- 
ory toward the low temperature dimerized antiferromag- 
netic transition. 

To conclude, below T* ~77 K, the inhomogeneous 
broadening of Cu NQR spectra and the impurity-induced 
Cu nuclear spin-lattice relaxation occur for undoped and 
more remarkably for Mg-doped CuGeOa. Precursory 
dimerization, inhomogeneous in real space, is suggested. 
The host antiferromagnetic correlation above Tgp is sup- 
pressed by Mg doping of a;=0.020. No magnetic critical 
divergence down to 4.2 K is a puzzle. 
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